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Abstract 

Electroweak onedoop radiative corrections to the process e + e~ — * ti are revisited. Two groups from 
Karlsruhe and Bielefeld/Zeuthen performed independent calculations of both (virtual and soft) QED 
contributions and weak virtual corrections. For the angular distribution an agreement of at least eight 
digits for the weak corrections and of at least seven digits for additional photonic corrections is established. 



* Work supported in part by the European Community's Human Potential Programme under contract HPRN-CT-2000- 
00149 Physics at Colliders. 

^E-mails: fl@physik.uni-bielefeld.de, hahn@feynarts.de, hollik@particle.uni-karlsruhe.de, Tord.Riemann@desy.de, 
cs@particle.uni-karlsruhe.de, Anja.Werthenbach@desy.de 



1 Introduction 



At TESLA, one of the most important production processes will be top-pair production well above the 
production threshold (in the continuum) 

e + (pi) + e~ (p 2 ) t (ps) + t (?4). (1) 

We expect several hundred thousand events, and the anticipated accuracy of the theoretical predictions 
to be a few per mille. What one has to calculate with such net precision is, of course, not only the 
two-fermion production process, but also the decay of the top quarks and the variety of quite different 
radiative corrections like real photonic bremsstrahlung, electroweak radiative corrections (EWRC), QCD 
corrections to the final state, and beamstrahlung. Potentially, new physics effects also have to be taken 
into account. 

The state of the art at the time of the public presentation of the TESLA project has been reviewed in 
Part III of the Technical Design Report Since then, the community has made progress in several 
directions. In this note we report on recent developments in the description of the electroweak radiative 
corrections to ([[]). Several studies on this topic are available in the literature, e.g. [||, |3|, ||, but no 
special effort was undertaken to perform a detailed numerical comparison between the results of different 
collaborations. Quite recently, a recalculation of the electroweak corrections has been performed by a 
Bielefeld/Zeuthen group ||, and this calculation was compared in detail with the Karlsruhe results.^] 

We summarize this comparison together with a few remarks on the underlying projects. 

2 One-loop EWRC to the reaction e + e~ — » ti 

In this section we introduce the one-loop differential cross-section of the process e + e~ — > ti in the 
Standard Model. 

To start with, we first present the differential Born cross-section, 

^Born /I »n2 



dcos# 



;i + cos 2 0)a£(s) + 2cosfl + (l-cos 2 #) 4' m (s), (2) 



with 



2 

4 ( s ) = ™ ^ [QlQl + 2Q e Q t v e v t Rexz(s) + {v 2 e + a 2 ) (v 2 + a*$) \xz(s)\ 2 ] , (3) 

2 

a FB( s ) = ~^T Ct $t 2QeQta e a t Rexz(s) + <iv e a e v t a t \xz(s)\ 2 , (4) 

2 

4' m (*) = [QlQ 2 t + 2Q e Q t v e v t Rexz(s) + {v 2 e + a 2 e )v 2 \ X z(s)\ 2 ] . (5) 



1 A series of numerical comparisons of the Bielefeld/Zeuthen group with the authors of JH) started in September 2001 and 
was later joined by the Karlsruhe group. 
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The Z propagator is contained in the factor 



XZ{S) ~ 4e 2 <fos 2 w s-m 2 z ' (6) 
ml = M 2 z -iM z F z , (7) 



and we use the following conventions 



Qe = "I, (8) 

a f = #(/) = ±\, (9) 
v f = a f (l-4\Q f \sm 2 8 w ), (10) 



together with 



ft = sli-*£. (") 



and the colour factor a = 3 of the top quark. Including the virtual electroweak corrections and the soft 
photonic contributions it proves to be convenient to use the notation of complex-valued form factors, 
each of them related to some helicity-based matrix element. For the production of massive fermion pairs, 
the cross-section depends on six independent form factors, 

^=W~\^\ 2 , (12) 

d COS 167T s 

with 



\M 1+Z \ 2 = {U 2 + T 2 ) (IF} 1 ) 2 + |Ff| 2 + |Ff | 2 + |Ff | 2 ) (13) 
+ 2 (U 2 - T 2 ) Re (F^Ff * + Ff F? 1 *) 
+ 2m 2 S (|Fl 1 | 2 -|Ff| 2 + |Ff 1 | 2 -|Ff| 2 ) 

+ 2^(UT-mt)Re(Fl 1 F\ 1 * + F| 1 Ff*) + 0{a 2 ) , 



with s = (pi + P2) 2 and 

T = S - (1 - A cos 6), (14) 
U = l{l + 0t cos 6). (15) 

The form factors appear in the decomposition of the matrix element according to 

M 1+z = F i ab M^ + Fl 1 M 3 ,u +Ff M 3M . (16) 

a,6=l,5 

The matrix elements A^3,n and ^3,51 emerge only in the calculation of the virtual corrections. Their 
interference with the Born amplitudes vanishes in the case of massless fermion pair production. The 
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other matrix elements M\^ a b receive additional contributions from soft photon corrections. In Born 
approximation, the non- vanishing form factors are 

F?b = (e 2 Xz(s) v e v t + e 2 Q e Q t ) / s, (17) 

F\% = -e 2 X z(s)v e a t /s, (18) 

Fl] B = -e 2 xz(s)v t a e /s, (19) 

Ff 5 B = e 2 X z{s)a e a t / s. (20) 

For completeness, we give also the matrix elements 

iMiab = [v(pi) -y^gaufa) ] [u(P3) In9bv(p4) ], a, b = 1, 5, (21) 

where g\ = 1 and g§ = 75. For massive loop corrections, two additional structures have to be added 

iM 3 ,u = [v(pi) tu(p 2 ) } Hps) tf 2 tv(p 4 ) ], (22) 
iM-3,51 = [v(pi) hlbu{p2) ] [Hps) ^-v{p 4 ) ]. (23) 

In Appendix A we give the collection of all Feynman diagrams relevant at one-loop level for the calculation 
in the 't Hooft-Feynman gauge. Ultra-violet (UV) divergences were treated within dimensional regulariza- 
tion. For renormalization the on-shell scheme was used 0. Infra-red (IR) divergences cancel analytically 
between the soft photon emission and the diagrams with virtual photons. This leads to a residual 
dependence of the cross-section on the maximal soft-photon energy, which has been arbitrarily chosen 
to be S™ ax = -^/s/10. Of course, after including real hard photonic bremsstrahlung this dependence 
disappears. 

In the on-shell renormalization scheme the weak mixing angle is defined by the weak boson masses. In 
this calculation the coupling constant c/ (|6|) is expressed by the electric charge and the weak mixing angle. 

M 2 

S m 2 6 w = (24) 
9 = -V- (25) 



2.1 The Karlsruhe Approach 

The Karlsruhe calculations have a long history dating back to the late 1980s, following the formulation 
worked out in |9|. We shall not describe them here in detail. They use the software packages FeynArts 
[10, 11], FormCalc [|l^], and LoopTools [12|, the latter of which is based on the FF library |l3|]. The 



calculations follow the scheme described in [14]: the diagrams are first generated with FeynArts, the 
resulting amplitudes are analytically simplified with FormCalc, whose output is then converted to a 
Fortran program which is linked with the LoopTools library to produce e.g. cross-sections. This process 
takes only a few minutes and is highly automated. 

Results of calculations performed in this way have been made available on The HEP Process Repository 
Web page |15[| . The actual codes used for comparison in this paper have been taken from there. In 
particular we used the files eett_sm.tar .gz and eett_smnoqed.tar .gz. 

We just mention that a file with the virtual one-loop corrections in the MSSM, eett_mssm.tar .gz, is 
also obtainable. 
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2.2 The Bielefeld/Zeuthen Approach 



The Bielefeld/Zeuthen group utilized the program package DIANA [16, 17|, which is a FORM interface for 
the package qgraph pl| for the creation of all the contributing Feynman diagrams to a given process. 



The output of DIANA is a sample of symbolic expressions which then were prepared for the integration 
of the loop momenta with FORM [19]. Two independent calculations which differ in many respects have 
been performed. In one of them, dimensional regularization is used for UV and IR divergences, while in 
the other the IR divergences are regularized by a finite, small photon mass. The resulting expressions in 
terms of Passarino-Veltman functions are calculated with a Fortran code, using two libraries of one-loop 
functions: the dimensionally-regularized version uses LoopTools [12] with minor additions, the other 
one the package FF JT^[. On-shell renormalization is performed following [20]. To eliminate the IR 
divergence of the virtual photonic contributions soft photon bremsstrahlung has been added; we followed 
the techniques described in [21]. 

The Fortran program topf it provides two possible outputs: Firstly the differential cross-section (and 
derived observables like total cross-section or forward-backward asymmetry), secondly six independent 
form factors in a given helicity basis, as introduced above, designed to be used in a Monte-Carlo program. 

It is worth mentioning that topf it is equipped to calculate the hard real photon corrections to the 
differential cross-section. 



3 Numerical Results 

The comparison of the two calculations is chosen to be performed on the level of the differential cross- 
section. We focus in particular on the following four contributions 



da 



d cos 



: Born cross-section which serves as a cross check of input parameters and conventions. 

Born 



da 
dcos 



: These numbers contain the Born cross-section plus the interference of Born (s- 

QED 

channel 7 and Z exchange) with one-loop virtual QED diagrams plus the absolute square of real 
soft photon radiation. No vacuum polarization diagrams, nor counter-diagrams, of the photon are 
taken into account, i.e. the running of the electromagnetic coupling constant is not included. It is 
assumed that this can be easily accounted for, including the variety of higher order corrections, by 
using the corresponding value of a. 



da 
dcos 



: The interference of Born (s-channel 7 and Z exchange) with one-loop virtual pure 

weak 

weak diagrams is shown. As before, the Born cross-section is added to the one-loop correction. The 
running of the electromagnetic coupling constant is included here. Renormalization is performed 
in the on-shell scheme and self-energy diagrams are taken into account as well. 



da 

dcos 



: All previous parts of the calculation are put together, i.e. the complete electroweak 

SM 



one-loop plus Born differential cross-section is given within the Standard Model. 



5 



All numbers in the following tables were obtained making use of the input parameters 



Tz = 


. 


a = 


9 

= 1/137.03599976, 


it 1 max 


rz /in 


M w = 


= 80.4514958 GeV , 


M z 


= 91. 1867 GeV , 


M H = 


-- 120 GeV 


m e = 


0.00051099907 GeV , 


m t 


= 173.8 GeV , 


m b = 


4.7 GeV , 


"V = 


0.105658389 GeV , 


m u 


= 0.062 GeV , 


m d = 


0.083 GeV 


m T = 


1.77705 GeV , 


m c 


= 1.5 GeV , 


m s = 


0.215 GeV 



(26) 



Effective quark masses reproducing the hadronic vacuum polarization contribution Aat with a sufficiently 
high accurancy have been chosen |22, [7]. As mentioned, the photon energy cut E™* is an arbitary 
quantity entering the final result if the radiation of hard photons in not taken into account. 

The numerical values of the differential cross-sections have been calculated at three typical Next-Linear- 
Collider energies. 

A /i = 500 GeV: 



COS I 



da 



d cos 9 



Born 



da 
dcos 9 



QED 



da 



d cos 



weak 



dcr 



dcos 9 



SM 



-0.9 



0.1088 3919 4075 
0.1088 3919 4075 



0.09866 4252 
0.09866 4253 



0.1242 59037 1 
0.1242 59037 6 



0.1140 8410 
0.1140 8410 



-0.5 



0.1422 
0.1422 



7506 9392 
7506 9392 



0.1285 0790 
0.1285 0790 



0.1568 
0.1568 



48371 9 
48371 8 



0.1430 8121 
0.1430 8121 



0. 



0.2254 
0.2254 



7046 4032 
7046 4033 



0.2023 9167 
0.2023 9167 



0.2402 
0.2402 



6680 4 
6680 3 



0.2171 8801 
0.2171 8801 



0.5 



0.3546 
0.3546 



6647 0332 
6647 0332 



0.3151 1724 
0.3151 1723 



0.3688 
0.3688 



8650 7 
8650 5 



0.3293 3727 
0.3293 3727 



0.9 



0.4911 
0.4911 



4371 5766 
4371 5767 



0.4307 1437 
0.4307 1437 



0.5033 
0.5033 



3751 2 
3750 8 



0.4429 0817 
0.4429 0816 



yfs = 700 GeV: 



cos 9 


da 


f da 


da } 


da 


LdcOS 6>J Bo rn 


Ldcos #J QED 


Ldcos 6M weak 


Ldcos 6>J S M 


-0.9 


0.05033 2867 3357 
0.05033 2867 3357 


0.04534 6950 
0.04534 6948 


0.05697 6919 7 
0.05697 6920 1 


0.05199 1003 
0.05199 1001 


-0.5 


0.06658 1166 2000 
0.06658 1166 2001 


0.05993 8980 
0.05993 8980 


0.07224 8491 7 
0.07224 8491 7 


0.06560 6305 
0.06560 6305 


0. 


0.1237 4052 8515 
0.1237 4052 8515 


0.1106 4932 
0.1106 4932 


0.1280 5635 7 
0.1280 5635 6 


0.1149 6514 
0.1149 6514 


0.5 


0.2218 4321 1646 
0.2218 4321 1646 


0.1955 2864 
0.1955 2864 


0.2224 6611 4 
0.2224 6611 3 


0.1961 5154 
0.1961 5154 


0.9 


0.3298 0454 9138 
0.3298 0454 9138 


0.2844 3984 
0.2844 3985 


0.3245 1514 8 
0.3245 1514 5 


0.2791 5044 
0.2791 5045 
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= 1000 GeV: 



cosO 


da 


da 


da 


da 


LdcOS 6 J Born 


Ldcos 6>J QED 


Ldcos 6>J weak 


Ldcos 0\ SM 


-0.9 


0.02278 5423 2732 
0.02278 5423 2732 


0.02036 5843 
0.02036 5844 


0.02552 1285 
0.02552 1286 


0.02310 1705 
0.02310 1706 


-0.5 


0.02978 2131 1031 
0.02978 2131 1031 


0.02674 1661 
0.02674 1663 


0.03186 3489 
0.03186 3490 


0.02882 3019 
0.02882 3021 


0. 


0.06118 0067 4224 
0.06118 0067 4224 


0.05453 9344 
0.05453 9344 


0.06159 1613 
0.06159 1613 


0.05495 0889 
0.05495 0889 


0.5 


0.1177 4694 9888 
0.1177 4694 9888 


0.1031 1627 
0.1031 1626 


0.1140 47686 
0.1140 47685 


0.09941 7009 
0.09941 6999 


0.9 


0.1811 2209 7086 
0.1811 2209 7086 


0.1540 3824 
0.1540 3823 


0.1713 46193 
0.1713 46191 


0.1442 6233 
0.1442 6232 



Table 1: The content of the tables is explained in the text, the upper and lower numbers correspond to 
the Karlsruhe and Bielefeld/Zeuthen approach, respectively. 

From Table 1 it is evident that the numerics of the weak virtual corrections are perfectly controlled with 
a gross agreement of the two approaches by at least eight digits. 

The pure photonic corrections agree to at least seven digits^]. The net result of the comparison being 
highly satisfactory, we did not push for more digits agreement. 

Finally, it has to be pointed out that we control here only what is often called the technical precision of a 
calculation, not to be confused with the precision of the prediction of some observable to be confronted 
with some realistic measurement. 



4 Summary 

We have demonstrated an agreement of seven to eight digits between the calculations of the Karlsruhe 
group and the Bielefeld/Zeuthen group for electroweak virtual and soft photonic corrections to the 
reaction e + e~ — > ti at the one-loop level. 

Naturally the claim of completeness of this calculation can only be applied to the electroweak corrections 
and only at one-loop accuracy. In particular the inclusion of QCD corrections and the finite life time 
of the top quarks have to be taken into account. To what extent higher-order weak corrections will be 
needed depends crucially on the expected experimental accuracy. Certainly, it is desirable to partially 
extend the calculations to higher-order effects, where a special emphasis could be given to large mass 



effects [0] and to the so-called Sudakov logarithms [24], which constitute the leading contributions at 
the next order in perturbation theory. Finally, needless to mention that various higher-order photonic 
corrections have to be added in high-precision numerical approach. 
2 The topfit part was analytically and numerically cross-checked against |23|. 
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Appendix 



A Feynman diagrams 

In cases where several particles are listed for a propagator, there is one diagram for each combination of 
particles. On fermion lines, an /j stands for {e, fx, r, u, c, t, d, s, b} and a for {v e , v^, v T }. 

Born diagrams: 




Counter-term diagrams: 




Self-energy diagrams: 

G w 
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Vertex diagrams: 




t t t t 



t t t t 




t t t t 
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